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ABSTRACT 

In this paper, we investigate the hnk between the hypervelocity stars (HVSs) discovered in the 
Galactic halo and the S-stars moving in the Galactic center (GC), under the hypothesis that they 
are both the products of the tidal breakup of the same population of stellar binaries by the central 
massive black hole (MBH). By adopting several hypothetical models for binaries to be injected into 
the vicinity of the MBH and doing numerical simulations, we realize the tidal breakup processes of 
the binaries and their follow-up evolution. Wc find that many statistical properties of the detected 
HVSs and S-stars can be reproduced under some binary injecting models, and their number ratio can 
be re-produced if the stellar initial mass function is top-heavy (e.g., with slope ~ —1.6). The total 
number of the captured companions is ~ 50 that have masses in the range ^ 3 — IMq and semimajor 
axes < 4000AU and survive to the present within their main-sequence lifetime. The innermost one 
is expected to have a semimajor axis ^ 300-1500AU and a pericenter distance ~ 10-200AU, with 
a significant probability of being closer to the MBH than S2. Future detection of such a closer star 
would offer an important test to general relativity. The majority of the surviving ejected companions 
of the S-stars are expected to be located at Galactocentric distances < 20 kpc, and have heliocentric 

radial velocities ^ — 500-1500kms^^ and proper motions up to ^ 5-20masyr^^. Future detection of 
these HVSs may provide evidence for the tidal-breakup formation mechanism of the S-stars. 
Subject headings: black hole physics — Galaxy: center — Galaxy: halo — Galaxy: kinematics and 
dynamics — Galaxy:structure 



1. INTRODUCTION 

More than a hundred young massive stars, mostly 
Wolf-Rayet/0 and B types, have been identified within 
a distance of ^ 0.5 pc from the massive black hole 
(MBH) in the Galactic Center (GCjiGillessen et al|[200l 
iLu et al.l [2009; ,Bartko et all l2010( ) . These young stars 
are empirically divided into the two groups: (1) the ma- 
jority of the young stars at a distance '-^ 0.04— (). 5 pc from 
the MBH are located on coherent disk-like structures, 
i.e., the clockwise rotating stellar (CWS) disk and the 
possible counter-clockwise rotatin g stellar (CCWS) disk 
(e.g.. iLevin fc Beloborodwl l2003t IPaumard et all [20061 : 
iLu et al.l 120091 ): and (2) the young stars within a dis- 
tance of 0.04 pc from the MBH (denoted as S-stars), ex- 
clusively B-dwarfs, are spatially isotropically distributed 
and their orbital eccentricities follow a distribution o f 
/e(e) (X e^-^ (e.g.. lGhez et al.ll2008l:[Giilessen et al.ll2009n . 
The existence of these young stars is quite puzzling as 
star formation in the vicinity of a MBH is thought to 
be strongly suppressed due to the tid al force from the 
MBH (i.e., the parad ox of youth, see iGhez etldl 120081 : 
iPaumard et al.|[2006[) . It is of great importance to ad- 
dress not only the formation of these stars but also the 
origin of their kinematics, which should encode fruitful 
information of the dynamical interplays between the cen- 
tral MBH and its environment. 

Young stars in the CWS (or CCWS) disk are proba- 
bly formed in a previously existing massive gaseous disk 
due to ins tabilities and fra.gmentation developed in it 
(•e.g..[Levin.2007,:.NaTOkshidl200a[Alexan et al.ll2008l: 
iBonnell fc Ricell2008| ). Young binary stars in the disk(s) 



may migrate or be scattered into the vicinity of the cen- 
tral MBH (e.g.. 'Madigan et al.ll2009D and then be tidallv 
broken up (e.g., .Hills, 1988; lYu fc Trema"hlil2003| ). One 
component of a broken- up binary may be ejected out as 
a hypervelocit y star (HVS) as discovered in the Galac- 
tic halo (e.g . [Brown et al.[ |2005l: lEdelmann et al.l [20051: 
iHirsch et al.l [2005[ ). and the other component may be 
captured onto a tig hter orbit similar t o tha t of the S- 
stars as proposed bv lGould fc QuillenI ([2003[ f^. If HVSs 
were initially originated from a stellar structure like the 
CWS disk, they may be spatially located c lose to the disk 
plane (|Lu et al.[[2010l: [Zhang et al.[[2010D . The current 
observations do show such a spatial correlation between 
the HVSs and the CWS disk, which suggests th at major- 
ity of the HVSs origina te from the CWS disk ([Lu et al.l 
[2?m[Zhm"g et al.l[20Tl . 

The HVSs discovered in the Galactic halo and the S- 
stars in the vicinity of the central MBH may naturally 
link to each other if they both are the products of the 
tidal breakup of stellar binaries in the vicinity of the 
central MBH. Therefore, it is interesting to simultane- 
ously investigate the properties of the HVSs in the Galac- 
tic halo (or the S-stars in the GC) and their captured 
(or ejected) companions, and probability distribution of 
these properties. Under the assumption that both the 
HVSs and S-stars arc the products of tidal breakup of 
stellar binaries, the working hypothesis in this paper, we 
construct a number of Monte-Carlo models to simulate 
the tidal breakup processes of stellar binaries in the GC 

^ For othe r me chanisms to produce the S-stars, see 
IBaruteau eTall l[20Tll) . 
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and check whether these models can accommodate the 
current observations, and make further predictions on 
both the companions of HVSs and that of S-stars for 
future observations Q 

This paper is organized as follows. In Section [21 we 
overview the tidal breakup processes of stellar binaries 
in the vicinity of a MBH and the dynamical connection 
between the ejected and captured components. Adopting 
relatively realistic initial conditions, we perform a large 
number of three-body experiments to realize the tidal 
breakup processes of stellar binaries in Section [3l As- 
suming a constant injection rate of stellar binaries into 
the vicinity of the central MBH and adopting the results 
from the three-body experiments on the ejected and cap- 
tured components, we use the Monte-Carlo simulations 
to produce both the HVSs and the S-stars. In Section HJ 
wc follow the orbital evolution of the captured stars to 
the present time by adop ting the autoregressi ve moving 
average (ARMA) model (jMadigan et al.ll2010[) . in which 
both the non-resonant relaxation (NR) and the resonant 
relaxation (RR) are included. The simulated S-stars ap- 
pear to be compatible with the observations of the S- 
stars. In Section [H we investigate the effects of different 
binary injection models on the number ratio of the sim- 
ulated HVSs to S-stars. The number ratio given by ob- 
servations can be reproduced if the initial mass function 
(IMF) of the primary components of stellar binaries is 
somewhat top-heavy. By calibrating the injection mod- 
els with observations, we estimate the number of the cap- 
tured (or ejected unbound) stars, as the companions of 
HVSs (or S-stars), that could be detected in the future. 
We also estimate the probability to have less massive 
stars captured on an orbit within that of S2 in Section [6l 
Conclusions are given in Section [T] 

For clarity, some notations of the variables that are 
frequently used in this paper are summarized in Table [TJ 
Given a physical variable X (e.g., mass, velocity, semi- 
major axis, eccentricity), the distribution function of X 
is denoted by fx{X) so that fx{X)dX represents the 
number of relevant objects with variable X being in the 
range X X + dX. 

2. OVERVIEW: TIDAL BREAKUP OF STELLAR BINARIES 
IN THE VICINITY OF A MBH 

A stellar binary may be broken up if it approaches a 
MBH within a distance of rtb = ab(3M,/m)^/'^, where 
Af, is the mass of the MBH, Ob is the semimajor axis of 
the binary, m = rUg + ml is the total mass of the binary, 
and nig and ml are the masses of the two components of 
the binary, respectively. During the breakup, one compo- 
nent of the binary, denoted as mg here, gains energy, and 
the other component ml loses energy. For an injecting 
stellar binary that is initially on a parabolic orbit relative 
to the MBH; the velocity of the binary mass center at its 

^ In principle, each HVS should have a companion left in the GC 
and each S-star should have a companion ejected to the Galactic 
halo. However, these companions could have left the main se- 
quence because of the limited lifetime and cannot be detected at 
the present time; and the captured companion of a HVS may even 
has been tidally disrupted by the central MBH and does not exist 
now. Considering of those cases, hereafter, the term 'companions' 
may have a broad meaning in that it includes the companions of 
previously existed HVSs or S-stars as well as those detectable at 
the present time; and the companions of HVSs and S-stars may 
have different numbers at the present time. 



periapsis to the MBH (~ rtb) is Vth ^ (GM,/rtb)^^^. 
The component mg receives a velocity change on the or- 
der of Svg ~ {ml /m) y^Gro/ob if the eccentricity of the 
stellar binary is 0, and it gains energy SE ~ mgVtb^Vg. 
The other component ml loses the same amount of en- 
ergy SE. If SE is sufSciently large, the component mg 
may manifest itself as a HVS with velocity at infinity 
Voo ^ y^2SE/mg if ignoring the deceleration due to the 
Galactic gravitational potential. The root mean square 
(rms) of Woo is approximately 



2,1/2 /0.1AU\ 
^ - Woo,0 ( j 



1/2 



2m.l 



1/2 



1/3 



M, 



4 X IO6M0 



1/6 



9iD), (1) 



where foo.o = 25 96km s ^ and g{D) is given by 
I Bromley et al.l (I2006D for injecting binaries on hyperbolic 
orbits with initial velocities at infinity of 250km s"-'^, i.e.. 



giD) =0.774 



0.02041? 
+7.62 X lO^^L*^ 
+8.62 X 10""L»^ 



6.23 X 10"*!?^ 

4.24 X 10"*L»'^ 



(2) 

where the penetration parameter D = lOOy-p ini/rtb char- 
acterizes the minimum distance where the binary ap- 
proaches the MBH; and rp.i„i is the initial pericenter 

1/2 

distance of the binary. The rms velocity (f^) appar- 
ently depends on the semimajor axis, the total mass and 
the mass ratio of the stellar binary, and the penetration 
parameter D. 

The exact value of Woo of the ejected component for any 
given stellar binary also depends on the relative orienta- 
tion of the stellar binary orbital plane to the orbital plane 
of the binary rotating around the MBH and the orbital 
phases of the two components at the time of its breakup. 

This dependence introduces a scatter of Voo around the 
1/2 

value (w^) given by Equation ([1]), as the orbital ori- 
entations of the injecting stellar binaries are probably 
random and the orbital phases of the two components 
are not fixed at the breakup time. Numerical simulations 
have shown that this scatter is approximately Gaussian 

with a dispersion of (t„^ ~ 0.2 (u^)^^^ (jBromlev et al.l 
I2006t iZhang et al.|[20l5l ). where the binary orbital ori- 
entations are assumed to be randomly distributed. The 
symmetry of the orbital phases of the two binary com- 
ponents (always at the opposite side to the mass center 
of the binary) ensures the same probability of receiv- 
ing energy for each star, which leads to the same ejec- 
tion probability for both components if the injecting bi- 
naries are initially on parabolic orbits (jSari "etall[2010l: 
iKobavashi et al.ll2012D . 

Stellar binaries on orbits bound to the MBH may ex- 
perience multiple close encounters with the MBH and 
the binary semimajor axes and eccentricities may be cu- 
mulatively excited to larger values until finally being 
broken up (jZhang et al.ll20l"0f ). The distribution of Voo 
for the ejected stars, produced during the first encoun- 
ters of the binaries with the MBH, follows a fitting for- 
mula similar to Equation 1^ over D ^ 20 — 150, i.e., 
g{D) oc 1 — (13/256)^, as the initial bounding energy of 
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TABLE 1 

Notation of some symbols 



Symbol 



Description 



M, 

rup 
rris 
m 
R 

^p.ini 

^oo ,ir 
Ani 

rth 

D 

a 

/3 

7 



mass of the central MBH 

mass of the primary component of an injecting stellar binary 
mass of the secondary component of an injecting stellar binary 
total mass of an injecting stellar binary, i.e., nip + rris 
nfis/mp 

initial scmimajor axis of an injecting stellar binary 

initial pericenter distance of the mass center of the injecting stellar binary to the MBH 
initial scmimajor axis of the orbit of an injecting stellar binary rotating around a central MBH 
initial velocity of the injecting stellar binary at infinity if the binary is on a hyperbolic orbit 
initial energy of the stellar binary 
tidal radius for the stellar binary 



orbital penetration parameter of the injecting stellar binary (e 
exponent of the power-law distribution of ab.ini 
exponent of the power-law distribution of ?'p,ini 
exponent of the power-law distribution of nip 



100rp,i„i/rtb) 



mass of the component that gains energy during the tidal breakup of a stellar binary 
mass of the component that loses energy during the tidal breakup of a stellar binary 
ml/nig^ 

exchange energy between the two components during the tidal breakup of a stellar binary 

mass of the ejected star alter the tidal breakup of a stellar binary 
mass of the captured star after the tidal breakup of a stellar binary 
velocity of the ejected component at infinity 
orbital scmimajor axis of the captured component 

orbital scmimajor axis of the captured component if the injecting binary is initially on a parabolic orbit 

orbital eccentricity of the captured component 

simulated total number of the ejected stars given a mass range 

simulated total number of the captured stars given a mass range 

simulated number of the detectable HVSs at the present time for given selection criteria 
simulated number of the detectable captured stars at the present time for given selection criteria 
simulated fraction of the ejected stars that survive to the present time on the main sequence 
simulated fraction of the captured stars that survive to the present time on the main sequence 
simulated fraction of the captured stars that have already been tidally disrupted until the present time 
simulated fraction of the captured stars that can be detected at the present time for given selection criteria 



ml 
1 

5E 



Hcap.O 
Ccap 

^^HVS 
cap 

luobs 
^ *cap 
pit 
-^HVS 
pit 

cap 
ptd 

cap 
E^obs 

cap 



the injectinK stellar bi naries is still significant (for details, 
see lZhang et al.ir201Q[ ). For multiple encounters, the en- 
ergy exchange 6E between the two components is deter- 
mined by the properties of the stellar binaries at the final 

1 /2 

revolutions. Our simulations show that (w^) of those 
ejected components for stellar binaries broken up within 
1000 revolutions around the MBH still follows Equation 
H]), but g{D) is now best fitted by 



,(i.)^l- 0.427( -I) -1.88(1^ 



-3.29 



D 
256 



1.37 



D 
256 



(3) 



for Z? < 300. 

For stellar binaries on bound orbits, the light compo- 
nent has a larger probability to escape away from the 
MBH because the specific energy it could gain is gener- 
ally larger than that of the heavy component in a coun- 
terpart case (see Eq. [ij. However, the difference in the 
ejection probability for the two components of the stel- 
lar binaries is significant only when the mass ratio of 
the massive ones to the light ones > 5 and 25E/mg is 
close to its initial bounding energy G'M,/2ab-»,ini, where 
ab-».ini denotes the initial scmimajor axis of the binary 
system compo sed of a steU ar binary and the MBH (see 
also iKobavaslii et al.l I2012D . 



The component ml of a broken-up stellar binary loses 
energy by an amount of SE and it is captured onto a 
tighter orbit with scmimajor axis Ccap- According to the 
energy conservation law, we roughly have 



1 2 



GmlM, 

2 (2 nan 



Si, 



(4) 



where Ei^i is the initial energy of the stellar binary, and 
it is ^ ^n^v"^ jjjj if the binary is initially on a hyperbolic 

orbit, or ^ ~GmM,/ (2ab-».ini) if on a bound orbit. The 
initial internal mechanical energy of the stellar binary 
— is ignored in Equation Qj. We now have the 

general form for Voo as v^o ^ ^2{SE + ^i?ini)/mg. If 
\Emi\ ^ 6E, the scmimajor axis of the captured star is 
GM, 



O'cap — '^cap,0 — Q 2 

nrr^r. ..J M, \ / 1 000km s" ^ \ ^ 

where q = ml/rrig. For the cases considered in this paper, 
the injecting stellar binaries are either initially on hyper- 
bolic orbits (but close to parabolic ones) or from stellar 
structures like the CWS disk, and thus Ocap ~ Ocap.o 
as approximately jE'inil ^ SE. Equation ^ shows the 
connection between the properties of the captured stars 
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left in the GC and that of their ejected companions in 
the Galactic bulge and halo. For those HVSs discovered 
in the Galactic halo with Voo ~ 700 — 1000km s~^I3 their 
companions left in the GC may be initially on orbits with 
semimajor axis in the range of ^ 3500 — 7000AU as stel- 
lar binaries with extreme mass ratios are rare. For the 
innermost S-star, i.e., the S2, of which the semimajor 
axis is ~ lOOOAU, its companion ejected out should have 
Woo ^ 1900kms"i if g - 1, and ~ eOOkms^^ if g - 0.1, 
respectively. For stellar binaries initially tightly bound 
to the MBH, we may have SE + 7^£^ini < 0, the compo- 
nent gaining energy either remains bound to the MBH 
or is ejected out with low velocities. 

The distribution of HVS properties is directly con- 
nected to the distribution of S-star properties (note that 
here we do not mean that an observed HVS in the Galac- 
tic halo is directly associated with an observed S-star in 
the GC as the products of the tidal breakup of the same 
binary star). According to Equation ([5]), the distribu- 
tion of the semimajor axis of the captured stars fa^^p 
is related to the distribution of the velocity of HVSs at 
infinity f^^ if m„ ^ ml, i.e.. 



fvaai'^Oo) '^Cap/ocap ("^Cap) 



_ GAf, 



(6) 



which suggests that any one of the two distributions 
above can be inferred from the other one. 

The velocity distribution of the ejected stars are mainly 
determined by the initial sets on the distributions of ab,ini 



and Tn 



since 



1/2 



'1/2 



(7) 



where g{D) denotes the dependence of the rms veloc- 
ity (w^)^^^ on the penetration parameter D, as shown 
in Equations ^ and for the cases of injecting bina- 
ries initially on hyperbolic orbits but close to parabolic 
ones and bound orbits like the stars in the CWS disk, 
respectively. The fitting forms of g{D) obtained from nu- 
merical experiments (Eqs. [2]and[3|) are decreasing func- 
tions in the range of 20 < D < 150 or 20 < £> < 300, 
and thus g{D) may be approximated as a monotonically 
decreasing function. We assume that the initial distri- 
bution of flbjiii and 

^p,ini are /'ab('^b,ini) ^ ini '^H*^ 

/rp(?'p.ini) fx fp respectively; and the probability of 
a stellar binary with semimajor axis ab,ini broken up by 
the central MBH at a penetrat ion distance -P is on ly a 
function of D, i.e., /d(-D) (see iBromlev et all [20061 ). If 

1/2 

q ^ \ and ignoring the scatter of Voo around (t^^) (i-e., 

1 /2 

^^oo ^ i^oo) ): the velocity distribution of the ejected 
components can be obtained as 



A„(Woo)0C 



d 



(X ■ 



/ob (Q-b,ini)/rp (''pani) X 
ini^^p.ini 



dVoo 

.ini 

d 



^ The estimated Voo for those detected HVSs depends on the 
Galactic potential model adopted, the val ues here are obt ained 
from the Galactic potential model given bv lXue et al.l II2008I '). 



CX V, 



(8) 

and this relation is valid only if g{D) is a monotonically 
decreasing function and it is independent of the detailed 
form of fD{D). Similarly, we also have 



lac 



, i"cap; ^ "cap ' 



(9) 



which is consistent with the simple relation given by 
Equation (jG]). The estimated slope of fv^{voo) (or 
/acap(^cap)) above is not affected by taking account of the 

1/2 

Gaussian- like scatter of Wqo around (w^) as the distri- 
bution is a power-law. If considering of the various mass 
ratios among the injecting binaries (see Eqs. [l]and[5]), 
however, the resulted slope of /^^ (vao) may be somewhat 
fiatter than the simple estimates above. Note also that a 
larger /3 may correspond to a slower migration/diffusion 
of stellar binaries into the low angular-momentum orbits 
or the vicinity of the central MBH, and lead to fewer 
HVSs at the high-velocity end and fewer captured stars 
in smaller distances to the MBH. 

The periapsis of a captured star ml is roughly ^ rtb 
and the characteristic eccentricity of the captured star is 



--cap 



1 



rth 



1 



^1/3(1 + ^)2/3 \^M, 



ml 



1/3 



(10) 



the Ccap depends on q and the mass ratio of the cap- 
tured star to the MBH. Considering of the Gaussian-like 
scatter in Voo and correspondingly the scatter in Ccap, 
the probability that the breakup of a stellar binary with 
given semimajor axis and mass of each component results 
in a captured star with eccentricity < ecap is roughly 



P{< ecap) 



1 



-erfc 



\J{1 - ecap)/(l - ecap) - 1 

V2a„^/(t;^)V2 



(11) 



where ecap is given by Equation ([TU]) and / {i 
0.2. To capture a n S2-like star (i.e ., Cc 



)l/2 ^ 

. , 0.887 and 

ml ^ 15Mq; see IGhez et"all 120081 and lCillessen et al.l 
I2009D directly through the tidal breakup of stellar bina- 
ries, it is necessary to have mg ^ ~ 15Af0 {q <C 1) 
and the probabilit y is ~ 0.25 if q = 0.1 according to 
Equation (HH) fsee iGould fc QuillenI [200l . The proba- 
bility to capture stars onto orbits with ecap < 0.8 is only 



10 



-17 



if m„ 



ml 



15Mq; and ~ 6 x lO"'^ even if 



10m/ = 15OAf0 (see also IGould fc QuillenI [2003). 
Since the number of S-stars in the GC is only on the order 
of a few tens, it is difficult to prod uce all the 9 observed 
S-stars with eccentricities < 0.8 (jGillessen et al.l [20091) 
directly by the tidal breakup of stellar binaries. In ad- 
dition, the captured stars may initially remain on a disk 
plane if their progenitor binaries are orig inated fro m disk- 
like stellar structure(s) as suggested bv ILu et "all (pOlOil . 
which is different from the isotopic distribution of the 
S-stars. Therefore, additional physical mechanism is re- 
quired to further make the captured stars evolve to or- 
bits with lower eccentricities and spatially isotropically 
distributed if the S-stars are really originated from the 
tidal breakup of binary stars. 

The re sonant relaxation processe s (RR), initially pro- 
posed bv lRauch fc Tremaind (jl996[ ). may cause the cap- 
tured stars dynamically evolving to their present or- 
bits as discussed by a number of authors fjLcvin. ,20071 : 
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iHopman fc Alexandeii 120061: iKoesis fc Tremaind 1201 ID . 
In Section|4j we will take into account the relaxation pro- 
cesses, including RR, to approximately follow the dynam- 
ical evolution of each "S-star" after its capture due to the 
tidal breakup of stellar binaries, we then check whether 
the eccentricity and spatial distributions of those surviv- 
ing "S-stars" are compatible with current observations. 
Note here we use the quotes around the term S-stars to 
represent all of those captured stars with mass in the 
range of ^ 7 — 15Mq; while the simulated S-stars (with- 
out quotes) represent those with mass ~ 7 — 15Mq sur- 
viving to the present time, which presumably correspond 
to the observed ones (see Section U]). 

3. MONTE-CARLO SIMULATIONS 

In this Section, we first adopt Monte-Carlo simulations 
to realize the tidal breakup processes and generate HVSs 
and "S-stars" , and then we investigate in detail the con- 
nection between the simulated HVSs and "S-stars" . We 
use the code D0RPI5 based o n the explicit fifth (fourth) - 
order Runge-Kut ta method ([Dormand fc Princi 119801 : 
iHaier et al.l 119931 ) to calculate the three-body interac- 
tions between a stellar binary and the c entral MBH. 
For de tails of the numerical calculations, see lZhang et alJ 
(I20T0I) . The successive dynamical evolution of the cap- 
tured stars in the GC and the kinematic motion of the 
produced HVSs in the Galactic potential will be dis- 
cussed in Sections 2] and m respectively. 

3.1. Initial Settings 

The mass of the central MBH is set to be 4 x IO^Mq 
through out the numerical calculatio ns in this paper 
(|Ghez et al.ll2008l: IGillessen et al.ll2009l) . 

For the injecting stellar binaries, the initial conditions 
are set as follows: 

• The distribution of the semimajor axes ab,ini 
follows the Opik-law, i.e., a = —1 (e.g., 
iKobulnickv fc FrveJ[2007l) . 

• The mass distribution of the primary stars mp fol- 
lows a power-law function, /mp(TOp) cx; m'^. The 
distribution of the secondary star (mg) or the mass 
ratio R = nis/m-p can be described by two pop- 
ulations: (a) a twin population, i.e., about 40% 
of binary stars have i? ~ 1, and (b) the rest bi- 
naries, w hich follow a dis tributio n of fR(R) ^ 
constant (IKobulnickv fc FryerM2007tlKiminki et al.l 
[2008l[2009li ! 

• The initial eccentricity of the injecting binary is as- 
sume d to be ejni = 0, as adopted in previo us works 
(e.g., iBromlev et al|[200l lAntonini et aI|[20ll FI. 

• The orientation of the inner binary orbital plane 
is chosen to be uniformly distributed in cos (j) for 
(/.G (0,^). 

For the orbits of the injecting stellar binaries relative to 
the central MBH, the initial conditions arc set as follows: 

Alternatively assuming the initial eccentricities ~ 0.3, the ve- 
locities of the resulted HVSs from the four models are roughly 
smaller than those obtained for ~ by < 10%. 



TABLE 2 

Different injection models for tidal break up of 
binaries 



Model 


7 


/3 


ab-.,mi{ Pc) 


foo,ini(kms 1)=^ 


Unbd-MSO 


-2.7 







250 


Disk-MSO 


-2.7 





0.04-0.5 




Disk-THO 


-0.45 





0.04-0.5 




Disk-TH2 


-0.45 


2 


0.04-0.5 





^ For the Unbd-MSO model, the injecting stellar binaries 
have initial velocities of 250km s~^ at infinity. 



• The stellar binaries are assumed to be initially in- 
jected from either disk-like stellar structures (simi- 
lar to the CWS disk) or infinity. If they were from 
structures like the CWS disk, the semimajor axes 
ab-.,ini follows a power-law distribution propor- 
tional to (X a^l'^i^-^ in the range of ~ 0.04 — 0.5 pc 
according to cu r rent observations on the CWS disk 
(iLu et al.|[200l iBartko et al.|[2009f ). If they were 
from infinity, i.e., unbound to the MBH, their 
initial velocities at infinity are set to Woo.ini = 
250kms-^ 

• If the injecting binaries were from disk-like stel- 
lar structures, the orientations of their orbits rel- 
ative to the MBH are assumed to satisfy a Gaus- 
sian distribution around the central planes of the 
stellar disks wi t h a standard deviat ion of 12° (cf., 
iLu et al.l [2OOI IBartko et all [20091) . The planes 
of the host disks are assumed to be the same 
as the two planes that best fit the observations, 
i.e., il,b) = (311°, -14°) and (176°, -53°), respec- 
tively, and these two planes are consistent with 
the CWS disk plane and the plane of the north- 
ern arm (Narm) of the mini-spiral in th e GC (or 
the outer warped part of the CWS disk; iLu et al.l 
120101 : iZhang et al.ll2010f FI. The injection rates from 
these two disks are assumed to be the same. 

• The periapsis that the injecting stellar binaries 
approach the MBH is simply assumed to follow 
a power law distribution, /rp(?'p,ini) oc r^i^i and 
/3 > 0. A larger value of (3 corresponds to a 
smaller fraction of the injecting stellar binaries that 
could approach the immediate vicinity of the cen- 
tral MBH. It is still not clear which mechanism is 
responsible for the migration (or diffusion) of stel- 
lar binaries into the vicinity of the central MBH, al- 
though the secular instability devel oped in a stellar 
disk is proposed to be a viable one (jMadigan et al.l 
l2009f ). Instead of incorporating the detailed mi- 
gration/diffusion process of the stellar binaries in 
the Monte-Carlo simulations below, we choose to 

^ Note that IBrown et al.l I I2012I ) recently reported five new un- 
bound HVSs discovered in the Galactic halo and re-analyzed the 
HVSs previously discovered. According to this new study, there 
are 17 unbound HVSs in the northern sky and the y are still consis- 
tent with being located on two planes revealed bv lLu et al.l II2010I ). 
That is, one of the disk plane is consistent with the CWS disk plane, 
while the other disk plane is more consistent with the warped outer 
part of the CWS disk and slightly deviates from the Narm plane. In 
this paper, we do not distinguish the Narm plane from the warped 
outer part of the CWS disk. 
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parametrize the migration/diffusion process quah- 
tatively by different values of /3 and a larger /? cor- 
responds to a slower migration/diffusion process. 

In this Section, we perform Monte-Carlo simulations 
by adopting four sets of initial conditions (as listed in 
Table [2]). In the first model, the stellar binaries are as- 
sumed to be injected from infinity with initial velocity of 
'^oo.ini = 250km s~^. For the primary components of the 
in.iecting bina ries, we adopt the Miller-Scalo IMF (e.g., 
iKroupal 120021 ). This model is denoted as "Unbd-MSO". 
For the other three models, the stellar binaries are as- 
sumed to be originated from stellar structures like the 
CWS disk, and the IMF of the primary components is 
either set to be the Miller-Scalo IMF or a top-heavy IMF 
with a slope of 7 = —0.45 as suggested by recent obser- 
vations of the disk stars fsee iBartko et al1l2010[ ). The 
slope of the initial distribution of the pericenter distance 
/3 is set to be either or 2. These models are denoted as 
"Disk-MSO" , "Disk-THO" and "Disk-TH2" , respectively 
The total number of three-body experiments is 10^ for 
each model with the initial settings described above. By 
comparing the results obtained from those different mod- 
els, one may be able to distinguish the effects of different 
settings on the IMF and the injection of stellar binaries. 

If not specified, those ejected or captured stars with 
mass in the range of ^ 3 — 15M0 are recorded, thus both 
the HVSs with mass ^ 3 — AMq and the captured stars 
with mass ~ 7 — 15Mq, corresponding to the currently 
detected ones, can be taken into account simultaneously. 
The ejected or captured stars with mass in the range 
of ~ 4 — 7Mq are also considered for completeness. For 
other ejected or captured stars with mass out of the range 
of 3 — I5M0, they may be either too faint to be detected 
or too massive with too short lifetime and thus with too 
small probability to survive. 

3.2. Numerical Results 

Figurc[T]shows the distribution of the tidally broken- up 
stellar binaries in the Woo — icap plane, where v^o is the ve- 
locity at infinity of the ejected component and Ocap is the 
semimajor axis of the captured component. As shown in 
panel (a), the majority of the simulated Voo — Ocap pairs 
obtained from the Unbd-MSO model are close to the one 
estimated from Equation ([5]) by setting ml /nig = 1 (solid 
line). The main reasons for this are: (1) the majority 
(70%) of the injecting stellar binaries have mass ratios 
q = Tnl/nig in the range of (1/2, 2) under the assump- 
tion of two populations set for the stellar binaries; and 
(2) all the injecting binaries have the same but negli- 
gible initial energy -Eini- A small number of v^o — Ocap 
pairs, which apparently deviate significantly away from 
the solid line (for q = 1 obtained from Eq. [5]) , are due 
to the breakup of the binaries with q substantially larger 
or smaller than 1 (below or above the solid line). For 
the other three models, the simulation results do not de- 
viate far away from the simple predictions by Equation 
(U) (for q ~ 1), except that fewer ejected stars at the 
high-velocity end are produced in the Disk-TH2 model 
than in the other models simply because not many stel- 
lar binaries can closely approach the MBH. The scatters 
of Voo around that predicted by Equation (g]) in pan- 
els (b)-(d) are more significant compared with that in 
panel (a), which is caused by one or the combination of 



the effects as follows: (1) a distribution of the negative 
initial energy of the injecting stellar binaries originated 
from stellar structure like the CWS disk (panels b-d); (2) 
relatively more progenitor binaries have q substantially 
larger or smaller than 1 in the cases with a top-heavy 
IMF (panels c and d); and (3) fewer stellar binaries ap- 
proach the immediate vicinity of the central MBH in the 
case of a large (3 (panel d). 

As seen from Figure [l] if the observed S-stars, with 
semimajor axes ^ 1000 — 4000AU, are produced by the 
tidal breakup of stellar binaries, their ejected compan- 
ions arc expected to have Voo ^ 1000 — 1600km s~^ in the 
Unbd-MSO model and 200 - 1500kms"i in the other 
models. The Disk-TH2 model produces fewer ejected 
stars with Voo substantially larger than 1000km s~^ com- 
pared with other models. The captured companions of 
the detected HVSs in the Galactic halo are more likely 
to have Ocap ~ 3000 - 8000AU in the Unbd-MSO model, 
which is consistent with the simple estimation by Equa- 
tion ©, and have Ocap ^ 1000 - 8000AU in the other 
models. The travel time of the detected HVSs from the 
GC to its current location is on the order of ^ 100 Myr, 
which suggests that their companions were captured 
^ 100 Myr ago and the orbits of the captured compan- 
ions may have been changed due to the dynamical inter- 
actions with its environment (see Section|3]). As shown in 
Figure [H for those captured stars with Ocap > lOOOOAU, 
the probability that they have ejected companions with 
Woo > 700 - 1000km s"^ is negligible. 

Figure [2] shows the distribution of those captured stars 
obtained from each model in the Ocap versus log(l — ecap) 
plane, where Ocap and ecap are their semimajor axes and 
eccentricities achieved right after they were captured, re- 
spectively. Relatively more captured stars with low ec- 
centricities are produced by the Disk-MSO model than 
by the Unbd-MSO model (see panels a and b) mainly 
because the injecting binaries can be broken up at rel- 
atively larger distance in the Disk-MSO model due to 
multiple encounters. And relatively more captured stars 
with low eccentricities and Ccap are produced in the Disk- 
TH model than that in the Unbd-MSO model (see panels 
a and c) because there are more injecting binaries with 
mass ratio q substantially less than 1 and ml ^ 7— 15Afo 
for a top-heavy IMF than that for the Miller-Scalo IMF. 
The Disk-TH2 model produces relatively more captured 
stars with smaller eccentricities for any given Ocap than 
the Disk-THO model (as shown in panels c and d), as 
those binaries are generally broken up at even larger dis- 
tances in the Disk-TH2 model because fewer binaries 
can approach the very inner region due to the steep- 
ness of the adopted /rp(''p,ini)- However, the eccen- 
tricities of those captured stars, even produced in the 
Disk-TH2 model, are still statistically significantly higher 
than that of the observed S-stars. The orbits of a num- 
ber of S-stars, including S2, can be directly produced 
in the Disk-THO model and the Disk-TH2 model if the 
injection rate of binaries is around a few times 10~^ to 
lO^'^yr"^ as set for those mod els (see similar rates ob- 
tained bv lBromlev et aLll2012D . According to Figure [U 
apparently it is extremely difficult to produce "S-stars" 
with Ccap < 0.8 directly through the tidal breakup mech- 
anism of stellar binaries in the vicinity of the MBH. Note 
also that fewer captured stars with Ccap < lOOOAU are 
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Fig. 1. — Number distributions of brokon-up stellar binaries in the Voc — Ucap plane, where v^c is the velocity at infinity of the ejected 
component and acap is the semimajor axes of its captured companion. Panels {a)-(d) show results obtained from the Unbd-MSO, Disk- 
MSO, Disk-THO, and Disk-TH2 models, respectively. The total number of the three-body experiments is 10"" for each model. The solid 
red line in each panel shows the estimation according to Equation Q for nig = ml by assuming I'oo ini = 250km s~^ in panel (a), 
and ai3_,,ini = 0.2 pc in panels (b)-(d), respectively. The dashed red lines above or below the solid lines represent the estimations for 
rrig/ml = 1/2 and 2, respectively. The magenta lines with triangles indicate the rms of Van for each bin of acap- The number of stars 
is counted in each of the acap and v^o bins (totally 25 X 25 bins with with bin size 400AU X lOOkms"^) and represented by the color 
brightness scales shown in the label. 

produced in the Disk-TH2 model compared with those 
in other models because stellar binaries are harder to ap- 
proach the innermost region than that in other models 
(see panel d in Figs. [1] and [5]). 

According to the simulations above, we find that ~ 
50% of the detected HVSs should have captured com- 
panions in the GC with mass ^ 3 — 4A'/q as shown in the 
left panel of Figure |31 And similarly ~ 60% of the ob- 
served S-stars should have ejected companions with mass 
'--^ 7 — 15Mq as shown in the right panel of Figure [3] To 
find the possible counterparts of those current observed 
S-stars and HVSs. we will focus on the ejected stars with 
mass 7~ 15Mq in the Galactic bulge and halo and the 
captured stars with mass ^ 3 — 4Mo in the GC. For com- 
pleteness, we also count the ejected and captured stars 
with mass ~ 4 — 7Mq produced in all the models (see 
Tab. HI). 

4. ORBITAL EVOLUTION OF THE CAPTURED STARS 

The orbits of the captured stars produced by the tidal 
breakup of stellar binaries may evolve due to dynami- 
cal interactions with the surrounding environments. In 
principle, two-body interactions between stars may cause 



exchanges of their angular momenta and energy. How- 
ever, the timcscale of the two-body relaxation 10^ yr) 
is too long for it to be effective in changing the or- 
bits of ca ptured stars within their main-sequence life - 
time (e.g., iHopman fc Alexandeii [200l lYu et al.|[2007[) . 
The resonant relaxation (RR) is an important dynami- 
cal process naturally resulted from the coherent torques 
between orbital averaged mass wires of stars moving in 
near-K eplerian potential proposed bv lRauch fc Trcmain^ 
([Ml), which can lead to changes in both the eccentric- 
ities (scalar RR) and orie ntations (vector RR) of star s 
moving m the GC (e.g.. iHopman fc Alexa^IdSl I2006D . 
The RR appears much more effective in changing the or- 
bital configuration of the "S -stars" than the non-resonant 
two-body r elaxation (NR: Hopman fc Alexander! 120061 : 
iPeretsI [2009: Kojisiii&LTremaine 2011^ Therefore, the 
scalar and vector RR may be crucial in the follow-up 
dynamical evolution of the orbits of the captured stars. 

The vector RR timescale is ~ 1 — 10 Myr in 
the region hosting the "S-stars", about one order 
of magnitude smaller than the scalar RR timescale 
(|Hopman fc Alexanded 120061 : lYu et al.l l2007| ) , and thus 
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Fig. 2. — Number distributions of the captured stars in the acap versus Iog(l — ecap) plane, where Ocap and ecap are the semimajor 
axes and eccentricities of the captured stars achieved right after their capture, respectively. Panels (a)-(d) show the results from the 
Unbd-MSO, Disk-MSO, Disk-THO, and Disk-TH2 model, respectively. The red open circles represent the observed S-stars with Ocap smaller 
than 4000 AU at the present time l IGillc sscn ct al. 200£i); and the magenta dashed and dotted lines are for the mean eccentricities given by 
Equation mOI I for (ml, q) = {IOMq , 1) and (IOMq , 1/10), respectively. The number of stars is counted in each of the acap and log(l — Ccap) 
bins (totally 25 X 25 bins with bin size 200AU X 0.08) and represented by the color brightness scales shown in the label. 




1 10 100 1 10 100 

mcap (Mo) (Mo) 

Fig. 3. — Mass probability distribution P(mcap) of the captured companions of the ~ 3 — AMq HVSs (left panel) and mass probability 
distribution P(moj) of the ejected companions of the ~ 7 — 15M0 captured stars (right panel). The solid (red), dotted (blue), dashed 
(magenta) and dot-dashed (cyan) lines represent the results obtained from the Unbd-MSO, the Disk-MSO, the Disk-THO, and the Disk-TH2 
models, respectively. 
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the captured stars can evolve to an isotropic distribu- 
tion on a timesca l e of ^ 10 Myr (iHopman fc Alexander! 
[200llPeretsl[2009HKocs"is fc Tremainell2011[ ) even if thev 
were originally on a plane-like structure. In this paper, 
we assume that the isotropic distribution of the S-stars 
can always be reproduced through the vector RR of the 
simulated captured stars within a timescale shorter than 
their lifetime. It has been suggested that the high ec- 
centric orbits of captured stars can dynamically evolve 
to that of the observed S-stars through the scalar RR 
within ~ 20 Myr (e.g., LPercts 2009). However, previous 
studies assume a simple distribution of the initial eccen- 
tricities and semima.jor axes of th e captured stars (e.g., 
iPeretd 120091 iMadigan et al.ll2010D . In this Section, we 
adopt the distribution of Ccap and Ocap resulted from the 
injection models studied in Scction|3l We follow the evo- 
lution of ecap and Ccap by taking both the RR and NR 
into account, and then compare the Ccap and Ocap dis- 
tributions of the captured stars surviving to the present 
time with that of the observed S-stars. 

We adopt the AR MA model first introduced by 
IMadigan et al.l (|2010( ) to perform Monte-Carlo simula- 
tions of the long-term evolution of the captured stars. In 
the ARMA model, the RR phase and the NR phase are 
unified, and the general relativistic precession of stars 
in the potential of the central MBH is also simultane- 
ously included. The ARMA model is characterized by 
the following three parameters: (1) the autoregressive 
parameter 0i; (2) the moving average parameter Oi; and 
(3) the parameter cti, which is the variance of a random 
variable e^^-* following the normal distribution. In the 
ARMA model, e^^-* represents the random walk motion 
of the NR phase. At a time step of one orbital period of 
a star, the variation in the absolute value of its angular 
momentum is 



and 



01 =exp(-— — ), 

t<t, = /0min[tpi.oc(a, e), iproc(a, e) 
1 



S 



1 + exp [-k[e- Ccrit)] 







ecrit(a, e) 



In A 



■ exp 




Ar 



M, 5tp 



iV<lnA 



Anr 

U = 0.52 + 0.62e - O.SGe^ + 0.21e^ - 0.29^^, 



(12) 

(13) 
(14) 
(15) 

(16) 



(17) 
(18) 

(19) 
(20) 



^ There is a mi sprint in the original form of S given by 
IJVLadigan et al.l 1 120101 . see their Eq. 33), i.e., the minus sign before 
k was missing. 



where 0.105, U = 1.2, fc = 30, ^Inr = 0.26, 

At — 1.57, A — M,/m,,, to* = lOM© is the averaged 
mass of the field stars, a and (5ip are the semimajor axis 
and orbital period of the star, e is the median value of the 
eccentricity of the field stars and it is \f\J^ for a thermal 
distribution, tp,.ec is the combined precession timescale 
for the Newtonian precession and the ge neral relativity 
prece ssion (see Eqs. 25, 27 and 28 in IMadigan et al.l 
I2010D . and A^< is the total number of stars within the 
radius equal to the semimajor axis of the captured star. 
Adopting a simple stellar cusp model, i.e., cx; 
we have iV< = A^h(''/''h)"^~", where rh represents the 
radius within which the mass of stars equals the MBH 
mass and A^h = A/,/m* = 4 x 10^ is the total num- 
ber of field stars within r\^. Similar to IMadigan et al.l 
(I2010D. we also assum e a Bahcall-Wolf cusp (a = 7/4, 
IBahcall fc WoHiri976f ) unless otherwise stated and corre- 
spondingly T\y = 2.3 pc. The variable superscript '(1)' 
in the above equations means that the time step is one 
orbital period of the star being investigated. In the time 
step of N period of the star [N ~ St/Stp), the model 
parameters {<j)N, Sn, cat) can be obtained from parame- 
ters of one period ((^i, 9i, ai). Fo r furth er details of the 
ARMA model, see IMadigan et al.l (|2010[) . 

The two-body non-resonant relaxatio n is also taken 
into a ccount in a way similar to that in IMadigan et al.l 
()2010[) . In a time step St, the energy change of a captured 
star due to the NR is given by 



AE = ^E 



St 

tNR 



1/2 



(21) 



where ^ is an independent normal random variable with 
zero mean and unit variance, <nr, is the NR timescale in 
the GC given by 



tNR = ^NR 



1 1 

InA" 



■St 



p- 



(22) 



Assuming a constant injection rate of stellar binaries 
over the past 250 Myr as suggested by the observations 
of the HVSs in the Galactic halo, we adopt the numerical 
results of the three-body experiments for each injection 
model in Section |3] and calculate the energy and angu- 
lar momentum evolution for each captured star by using 
the ARMA model. In these calculations, we also take 
account of the effects of the limited lifetime of the cap- 
tured stars on the main sequence and the tidal disruption 
of those captured stars moving too close to the MBH. We 
remove those captured stars if they move away from the 
main sequence or approach the MBH within a distance of 
rp < r"i, where r^^ = {2M,/mY/^R^ and is the stel- 
lar radius. Finally, we obtain the present-day semimajor 
axis and eccentricity distributions of the captured stars, 
which can be used to be compared to the observational 
distributions and constrain the models. 

The left panel of Figure |4] shows the cumulative dis- 
tributions of flcap of the captured stars surviving to the 
present time and that of the observed S-stars. The thick 
lines represent the captured stars with mass ^7 — 15Mq 
surviving to the present time (the simulated S-stars), 
roughly corresponding to the observed S-stars in the GC, 
and the thin lines represent the captured stars with mass 
~ 3 — 4M0, roughly corresponding to the captured com- 
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Fig. 4. — The cumulative distributions of the semimajor axis (left panel) and eccentrici ty (right panel) of th e captured stars surviving 
to the present time. The histograms represent the distributions of the observed S-stars IIGillessen et aLIPOOgD . The solid (red), dotted 
(blue), dashed (magenta) and dot-dashed (cyan) curves show the results obtained from the Unbd-MSO, Disk-MSO, Disk-THO, and E)isk-TH2 
models, respectively. The thick and thin curves represent the distributions of those captured stars with mass 7 — 15Mq and 3 — 4Mq, 
respectively. In the right panel, the dot-dot-dashed line represents the cumulative eccentricity distribution proportional to ej?^p as suggested 
by the observations IIGhez et al.|[200g : [Hillessen et al.ll2009l ). 

panions of those HVSs detected in the Galactic halo. Al- 
though the lifetime of less massive stars on the main 
sequence is substantially longer and thus the dynamical 
evolution time is longer than that of the massive ones, 
the cumulative distribution of acap of the light captured 
stars is only slightly different from that of the massive 



thus longer dynamical evolution time, and Equation ([5]) 
is no longer reliable to provide estimations on the velocity 
of the ejected companions of those less massive captured 
stars by using their present-day semimajor axes. 

The right panel of Figurc|4]shows the cumulative eccen- 
tricity distributions obtained from different models. As 
ones (see left panel of Fig. |4|). The slope of the Ccap dis- seen from Figure[4l the Ccap distributions are only slightly 



tribution is affected most by the /rp(''p,ini) distribution. 
Relatively fewer captured stars are produced in the in- 
ner region by the Disk-TH2 model compared with that 
obtained by the other models, and the fraction of cap- 
tured stars with semimajor axes < Ccap is proportional 



to 



,2.0 



for the Disk-TH2 model but to 



,10-1. 5 



for 



cap ^ v-^iv.. ^^^^^ ^ ^^^^^^^ "'cap 

the other three models. For those models with /3 = 0, 
the flcap distributions obtained from the numerical simu- 
lations is roughly consistent with the simple estimations 
from Equation ([9]) (e.g., the slope is 1 for a = —1 and 
f3 = according to Eq.[S]). For the Disk-TH2 model, how- 
ever, the Ccap distribution seems flatter than the simple 
expectation, i.e., a slope of 3 (for a = — 1 and /3 = 2). 
The flatter slope of /(ocap) resulted from the Disk-TH2 
model may be due to the effect of various mass ratios of 
the injecting binaries, which is included in the numeri- 
cal simulations but ignored in deriving Equation ^ (see 
panel d of Figure [1]). The Kolmogorov-Simirnov (K-S) 
tests find the likelihoods of 0.02, 0.2, 0.05, and 0.15 that 
the flcap distribution of the observed S-stars is the same 
as that of the simulated S-stars for the four models, re- 
spectively, which suggests that the latter three models 
may be compatible with the observational acap distribu- 
tion. 

For the majority of the simulated S-stars (^ 7— I5M0), 
the relative changes in their energy due to dynamical re- 
laxation after their capture are less than ^30% and their 
semimajor axes do not deviate much from the initial val- 
ues right after their capture. Therefore, the distribution 
of the semimajor axis of currently observed S-stars can 
provide some information on their ejected companions 
(see Eqs. |4]and[5]). For those captured stars with mass 
~ 3 — 4M0, however, their relative energy changes can be 
as large as 1 mainly because of their longer lifetime and 



different for different injection models because the initial 
Ccap are all close to 1 in all the models. For those cap- 
tured stars with mass ~ 3— 4M0, their present eccentrici- 
ties are relatively lower than that of the observed S-stars 
because of their longer dynamical evolution time. For 
those simulated S-stars (with mass ^ 7 — I5M0) resulted 
from any of the four injection models, their ecap distribu- 
tion is similar to that of the observed S-stars. The K-S 
tests find a likelihood of ^ 0.1 — 0.6 that the eccentricity 
distribution of the observed S-stars is the same as that of 
the simulated S-stars for all the four models. However, 
it appears that the simulations slightly over-produce the 
stars with high eccentricities (ccap close to 1) with respect 
to the observations because of the ineffectiveness of the 
RR process for th ose captured star s with extremely high 
eccentricities (Ma digan et al.l[2010( ) . One possible expla- 
nation to this inconsistency between the observations and 
simulations could be the observational bias in detecting 
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ine timescaie 01 tne fiK p rocess also aepends 
on the mass of the fie ld stars (|Madigan et al.l 120101 : 
iRauch fc Tremailii ll996'). To check this dependence, we 
perform additional simulations by setting the mass of 
field stars to m^, = 5AIq or = 2OM0 but with the to- 
tal mass of the field stars fixed. According to the results 
of these simulations, we find that the simulated S-stars 
are on orbits with too high eccentricities compared with 
the observational ones if = 5Mq field stars, or on 
orbits with too low eccentricities if m^, = 20Mq. 

The timescaie of the RR process becomes much longer 
at the distance of the disks. Stars in this region are less 
affected by the relaxation processes and may well pre- 
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serve some of their initial orbital configurations. Obser- 
vations find many B-dwarfs in disk regions, with high ec- 
centricities and more extended spatial distribution than 
disk stars. The resulted eccentricity-distance distribu- 
tion of these stars by the RR compared with the obser- 
vations may p rovide useful const raints on the formation 
of the S-stars (jPerets et al.l 120101 ). Our simulations also 
produce many B-dwarfs in the disk region and their ra- 
dial distribution is similar to the initial input ones for 
the injecting binaries. However, there should also exist 
B-dwarfs in the disk region that are initially formed as 
single stars and the fraction of these stars is not clear 
yet. A detailed dynamical study for the B-dwarfs in the 
disk region is complicated and beyond the scope of this 
paper. 

5. EJECTED HVSS IN THE GALACTIC BULGE AND HALO 

The ejected stars move away from the GC after the 
breakup of their progenitor binaries and their velocities 
are gradually decelerated in the Galactic gravitational 
potential. Some of them are unbound to the Galactic 
potential and can travel to the Galactic halo and may 
appear as the detected HVSs if their main-sequence life- 
time is long enough compared to the travel time; while 
others with lower ejecting velocity may return to the GC. 
To follow the subsequent motion of the ejected compo- 
ne nts, we adopt the Milky Way potential model given 
bv IXue et al.l (|20Q8| ). which involves four components, 
including the contributions from the central MBH, the 
Galactic bulge, the Galactic disc, and the Galactic halo, 
i.e., 

$ $BH + $bulgc + $disk + $halo, (23) 



where 



$BH = -GA/./r, 



Igc 



Igc 



7-buli 



disk 



halo 



ln(l 



-), 



(24) 
(25) 



(26) 



(27) 



respectively. The model parameters for the last three 
components are Mbuigc — 1-5 x IO^^Mq, Afdisk = 5 x 
10^°Afo, the core radius rbuigc = 0.6 kpc, the scale length 

6 4 kpc, and ps = | i^^l^cyZfii+c) ' where Pc is the cos- 
mic critical density, Avir = 200, fim is the cosmic fraction 
of matter, the virial radius r^ir = 267 kpc, and the con- 
centration c = 12. The bulge, the disk, and the halo po- 
tentials adopted here are all spherical. If adopting non- 
spherical potentials, i.e., a triaxial bulge/halo and a flat- 
tened disc potential, the bending effect due to the non- 
spherical component on the trajectories of ejected stars 
is important only for those with Voo ^ 400km s~^ on a 
timcscale of > 500 Myr, but it is negligible for HVSs with 
relatively high speeds (see lYu fc Madaul I2007D . Note 
that the radial distribution of the ejected stars surviving 
to the present time (and correspondingly the predicted 
number of the detectable HVSs) may be slightly different 
if adopting a different Galactic potential model. 



The total number of detectable HVSs depends directly 
not only on how many stellar binaries can be injected 
into the immediate vicinity of the MBH, but also on 
the lifetime of these stars and the detailed settings on 
the IMF, semimajor axis, and periapsis of the inject- 
ing stellar binaries (see Section [3]). If the stellar bi- 
naries are injected from a far away region, the injec- 
tion rate can be est i mated through the loss cone the- 
ory (|Yu fc Tremaind [2001 iPeretsI 120091 ): if the inject- 
ing ste llar binaries originated from central stellar disks 
(jLu et al. 2010; Zhang et al. 201(|), the injection rate is 
difficult to estimate as the mecha nism responsible for i t 
is still not clearly understood (cf. iMadigan et al.|[20Q9D . 
In principle, it is plausible to observationally calibrate 
the injection rate by the numbers of the detected HVSs 
and S-stars. But this calibration becomes complicated 
if considering of the uncertainties in the settings of the 
distributions /ab(ab,ini), /rp(''p,ini), and IMF of the in- 
jecting binaries, etc. 

5.1. The numbers of the HVSs/S-stars and their 
number ratio 

The 3 — 4Af0 HVSs detected in the Galactic halo and 
the S-stars in the GC should be linked to each other 
under the working hypothesis of this paper. The total 
numbers of the simulated 3 — 4Af0 HVSs and 7 — 15Mq 
S-stars depend not only on the injection rate of stellar 
binaries but also on the detailed settings on the injec- 
tion models. However, the number ratio of the simu- 
lated 3 - 4A'/0 HVSs to the simulated 7 - 15Af0 S-stars 
may depend only on the details of the injection models 
described in Section [31 but not on the injection rate of 
stellar binaries. Any viable model should produce a num- 
ber ratio compatible with the observations on the HVSs 
and the S-stars, and thus this number ratio may provide 
important constraints on the models. 

We obtain both the numbers of the simulated HVSs 
and the captured stars surviving to the present time and 
their number ratio by Monte-Carlo simulations as fol- 
lows. First, we obtain the total number of initially cap- 
tured (or ejected) stars, N^'^^ (or A^hvs)' with mass in 
the ranges of 3 — 4Af0, 4 — 7Af0, and 7 — 15Af0, respec- 
tively. To do this, we assume a constant injection rate 
of binaries and randomly set the injection events over 
the past 250 Myr, and for each injection event we ran- 
domly assign it to a three-body experiment conducted 
in Section [3| and adopt the results from the experiment. 
Second, we consider the limited lifetime of each ejected 
and captured star, the motion of each ejected star in the 
Galactic potential, and the dynamical evolution of each 
captured star, and then obtain the fraction of the ejected 
stars (^Hvs) that still remain on the main sequence at 
the present time or the similar fraction for the captured 
stars (-Fcap)j and the fraction of the captured stars (i^cap) 
that have already been tidally disrupted until the present 
time. Third, we consider the kinematic selection criteria 
and obtain the fractions of those ejected and captured 
stars according to given observational selection criteria, 
i.e., i^HVS ^^'^ ^cap; respectively. The selection crite- 
ria are similar to those adopted in selecting the detected 
HVSs and S-stars, i.e., an ejected star is labeled as a 
detectable HVS if its heliocentric radial velocity in the 
Galactic rest frame is |wrf | > 275kms~^ or its proper mo- 
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TABLE 3 

The number ratio of the simulated 3 — 4Mq HVSs to the simulated 

S-STARS. 



Model 


7 


P 


jytot 


pit 

cap 


ptd 

cap 


T-iobs,rf 

-'^hvs 


^cap 


oba.rf 

"hvs 


Unbd-MSO 


-2.7 





3.0 


8.8 


0.59 


0.27 


0.51 


27 


Disk-MSO 


-2.7 





1.8 


8.9 


0.51 


0.15 


0.41 


12 


Disk-THO 


-0.45 





0.22 


10 


0.47 


0.19 


0.54 


1.4 


Disk-TH2 


-0.45 


2 


0.27 


9.9 


0.34 


0.06 


0.21 


1.3 


Disk-IMO 


-1.6 





0.61 


9.5 


0.50 


0.16 


0.47 


4.0 


Disk-IM2 


-1.6 


2 


0.75 


9.4 


0.22 


0.01 


0.03 


2.8 



Note. — The N^^^ and Af*"* represent the total number of the ejected 
stars with mass 3 — 4M0 and the total number of the captured stars with 
mass 7— 15Mq that are generated by the tidal breakup of stellar binaries in 
the GC for each model, respectively; the F^^^ and -F^ap denote the fraction 
of the ejected and captured stars that still remain on the main sequence of 
their stellar evolution at the end of our simulations; the -F'^p denotes the 
fraction of the captured stars that have been tidally disrupted by the central 

MBH before the end of our simulations; the F^yg^' denotes the fraction of 
the ejected stars appear as the detected HVSs, where an ejected star is 
taken as a detectable HVS if its heliocentric radial velocity in the Galactic 
rest frame is \v-^f\ > 275kms~-', its velocity at infinity is Vaa > 750kms~-'^, 
and its distance from the GC is in the range of 40 — 130 kpc; /°^p denotes 
the fraction of the captured stars that are within radii < 4000AU from the 

MBH; and the N"^'^^ is the total number of those detectable HVSs with 
mass 3 — A.Mq in the Galactic halo and N°^^ is the simulated number of 
the captured stars with mass 7 — 15Mq in the GC, which correspond to the 
observed ones. 



TABLE 4 

Predicted numbers in different models 



Model 


7 


P 


injection rate 
(10-^yr-l) 




3 - 4Mq 






4 - 7Mq 






7 - 15AfQ 




A TObs,rt 


II T-obSjpm 


M-obs 
cap 


11 robs.rt 
^^HVS 


11 T-obSjpm 


yyobs 
' cap 


II rObs,rt 

^^hvs 


11 i-obSjpm 


^bbs 
cap 


Unbd-MSO 


-2.7 





13 (2.2) 


504 (79) 


189 (30) 


86 (13) 


178 (28) 


58 (9) 


33 (5) 


60 (9) 


33 (5) 


17 (3) 


Disk-MSO 


-2.7 





7.0 (2.7) 


177 (79) 


66 (29) 


82 (36) 


62 (28) 


21 (9) 


32 (14) 


29 (13) 


17 (8) 


17 (8) 


Disk-THO 


-0.45 





1.8 (5.8) 


23 (79) 


8 (27) 


11 (37) 


28 (97) 


8 (30) 


11 (39) 


29 (103) 


19 (66) 


17 (60) 


Disk-TH2 


-0.45 


2 


6.2 (23) 


18 (79) 


7 (31) 


8 (34) 


16 (71) 


5 (23) 


9 (40) 


20 (87) 


15 (64) 


17 (75) 


Disk-IMO 


-1.6 





3.1 (3.6) 


61 (79) 


22 (29) 


28 (36) 


37 (48) 


12 (15) 


18 (23) 


29 (38) 


18 (23) 


17 (22) 


Disk-IM2 


-1.6 


2 


60 (100) 


48 (79) 


21 (35) 


32 (52) 


24 (39) 


7(11) 


18 (29) 


20 (33) 


14 (23) 


17 (28) 



Note. — The numbers of the simulated detectable HVSs and the captured stars surviving in the GC at the present time in different mass 
ranges, obtained from different models. The injection rate of stellar binaries is assumed to be a constant over the past 250 Myr, which enables 
the production of 17 simulated S-stars [or 79 unbound 3 ~ 4Mq HVSs (numbers in the brackets)]. The Af^yg'^ denotes the number of HVSs with 
\vj-f\ > 275kms^^ and Vaa > 750kms^^; and the A^yg''™ denotes the number of the HVSs with proper motion > 5masyr~^ in the heliocentric 
rest frame and Vao > 750kms~^. In order to compare to the observations, a simulated HVS with mass of 3 — 4Mq is counted as a detectable 
HVS additionally if its distance to the GC is in the range of 40 — 130 kpc. The A°^p is the total number of the captmed stars surviving in the 
GC with present-day flcap < 4000AU in the simulations. 

stars and consequently the number ratio of these two 
populations. (1) The difference between the lifetime of 
the simulated HVSs and S-stars: the detected HVSs are 
in the mass range of ~ 3 — 4Mq, which are substantially 
smaller than that of the observed S-stars (~ 7 — ISMq) 
in the GC. The lifetime difference leads to an enhance- 
ment in the number ratio of the detectable HVSs to the 
S-stars roughly by a factor of 10 under the assumption 
of a constant rate of injecting stellar binaries into the 
vicinity of the MBH over the past 250 Myr. (2) The 
place where the stellar binaries are originated: in the 
"Unbd-MSO" model, relatively more stars with high ve- 
locities (e.g., > 1000km s^^) are generated than those 
in the other models. (3) The distribution of pericenter 
distance of those injecting stellar binaries, which is re- 
lated to the speed of the migration or diffusion process 



tion in the heliocentric rest frame is > 5masyr~^, and 
its velocity at infinity is Voo > 750km s^^, and captured 
stars with semimajor axis < 4000AU are counted as de- 
tectable S-stars. For those simulated HVSs with mass 
3 — 4M0, we put an additional cut on their distances 
from the GC, i.e., from 40 to 130 kpc, in order to com- 
pare them to current observations. Finally, we obtain the 
number ratio of the detectable HVSs to the detectable S- 
stars for each model as 

/\7obs Ar^ot pit pobs 

^^HVS _ ^^HVS ^ -^HVS ^ -^HVS f2a\ 

KTohs /Vtot pit (1 — Ptd Pobs ' ^ ' 

cap cap cap \ cap/ cap 

The simulation results are listed in Tables |3] and |3] for 
each model. 

Here we comment on a few factors that affect the pre- 
dicted numbers of the simulated detectable HVSs and S- 



S-stars in the GC and HVSs in the Galactic halo 
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of those binaries to the immediate vicinity of the central 
MBH: a change in this distribution may result in either 
a significant increase or decrease in both the number of 
HVSs and that of captured stars, but their number ratio 
is not affected much. 

The complete survey conducted by iBrown et alJ 
(|2009a[) has detected 14 unbound HVSs with mass ~ 
3—AMq in a sky area of ^ 7300 deg^ and the total num- 
ber of similar HVSs in the whole sky should be '-^ 79±2lQ 
Observations have revealed 17 S-stars within a distance 
of ~ 4000AU from the central MBH. The number ra- 
tio of the detected HVSs to the S-stars is ~ 3.4 — 5.9. 
The number ratio resulted from any of the top four in- 
jection models listed in Table [2] is inconsistent with the 
observational ones. Both the Unbd-MSO model and the 
Disk-MSO model give a number ratio substantially larger 
than that inferred from observations, while the Disk-THO 
model and the Disk-TH2 model give too small number 
r atios. 

iFiger et all (jl999D suggest that the IMF of young star 
clusters in the GC, i.e., the Arches cluster, may be top- 
heavy and the IMF slope is 7 ^ —1.6, although some 
later studies argued that the IMF of Arches cluster may 
be still consistent with the Sal peter one by assumiii g 
continuous star formation (e.g.. iLockmann et al.|[20Tol ). 
Considering of this, we adopt an IMF with a slope of 
7 = —1.6 and perform two additional injecting models, 
i.e., "Disk-IMG" and "Disk-IM2", as listed in Tables [3] 
and|4l Our calculations show that the number ratios pro- 
duced by the two models are close to the observational 
ones. Obviously, the number ratio of the simulated HVSs 
to S-stars is significantly affected by the adopted IMF. 
Adopting a steeper IMF may lead to a larger number 
ratio. 

The velocity distribution of the detected HVSs sug- 
gests a slow migration/diffusion of stellar binaries into 
the immediat e vicinity of the cen tral MBH (see detailed 
discussions in lZhang et al.ll2010D . The Disk-IM2 model 
can produce a velocity distribution similar to the ob- 
servational ones, and the large /3 adopted in this model 
also suggests a slow migration/diffusion of binaries into 
the vicinity of the central MBH. However, all the other 
models appear to generate too many HVSs at the high- 
velocity end and thus a too flat velocity distribution com- 
pared with the observational ones. This inconsistency 
could be due to many factors. For example, (1) the ve- 
locity distribution of the detected HVSs could be biased 
due to either the small number statistics and/or the un- 
certainties in estimating the observational selection ef- 
fects; (2) the uncertainties in the initial settings of the 
injecting binaries could also lead to some change in the 
velocity distribution, e.g., the velocity distribution may 
be steeper if the semimajor axis distribution of stellar 
binaries is log-normal, rather than follow the Opik-law 

One sdO type star with mass ~ IMq and another massive 
HVS in the southern hemisphere with mass ~ 9Mq in the survey 
are not included in the number. 

^ Considering of th e new results on searching HVSs reported 
bv lBrown et ahf | |2012| ). this number could be slightly higher, i.e., 
95 ± 23. And if assuming that the detected HVSs were originated 
from two disk-like stellar structure s, i.e., the CWS disk plane and 
the Narm plane, as suggested bv ILu et al.l II2010I ). the expected 
total number of HVSs with mass ~ 3 — 4Mq is 75 ±28. The number 
ratios of the detected HVSs to the S-stars become ~ 2.2 — 5.0. 



(jSesana et al.l 120071: [Zhang et al.ll2010D : and (3) the in- 
jection of binaries to the vicinity of the central MBH 
may be quite different from the simple models adopted 
in this paper because of the complicated environment of 
the very central region, e.g., the possible existence of a 
number of stellar-mass BHs or an intermediate-mass BH 
within lOOAU. 

For each model, the injection rate of stellar binaries 
can be calibrated to produce the numbers of the ob- 
served HVSs and S-stars. Under the initial settings for 
each model described in Section [3l the numbers of the 
simulated HVSs and the S-stars, similar to the detected 
ones, are listed in Table H] The calibrated injection rates 
are also listed in Table |4] for each model and they are 
roughly on the order of ~ 10~^ — 10~^ yr~^. However, 
one should be cautious about that this injection rate de- 
pends on the initial settings on the distributions of the 
semimajor axes and the periapses of the injecting bina- 
ries. According to our simulations, many of the injected 
binaries are not tidally broken up, and the breakup rate 
is a factor of 3 — 6 times smaller than the injection rate 
for those models studied in this paper, i.e., about a few 
times 10~^ — 10~^ yr~^, which is consistent with the es- 
timates bv lBromley et al1l2012l 

5.2. The ejected companions of the S-stars 

The HVSs discovered in the Galactic halo are typically 
in the mass range ^ 3 — AMq . HVSs with other masses 
should also be populated in the Galactic bulge and halo. 
In this Section, we investigate the properties of the simu- 
lated HVSs with mass ~ 7 — ISM©, which are most likely 
to be the companions of the "S-stars" . 

Panel (a) of Figure [5] shows the cumulative distribu- 
tion of the Galactocentric distances of the high-mass 
(^ 7 — 15M0) HVSs. The predicted numbers of these 
HVSs from different models are listed in Table |31 The 
total number of the simulated unbound HVSs, as com- 
panions of the population of the S-stars in the GC, is 
~ 20 — 30 according to the Disk-IM2 model, which is 
able to reproduce the numbers of the observed HVSs 
and S-stars. The majority of these unbound HVSs are 
at distances of a few to a few tens kpc from the GC, 
which are much closer to the GC than the detected 
3 — 4A/q HVSs mainly because a high-mass star has a 
shorter main-sequence lifetime and thus the distance it 
can travel within the lifetime is small. The close dis- 
tances of these HVSs from the GC suggest that the veloc- 
ity vectors of many of them are not along our line of sight. 
Their three-dimensional velocities range from 500km s~^ 
to 2000km s~^ in the Galactocentric rest frame (see panel 
(b) in Fig. [5]); and their heliocentric radial velocities 
in the Galactic rest frame range from — 500km s~^ to 
1500km s~^ (see panel (c) in Figure [5]), where the nega- 
tive and positive velocities represent moving toward and 
away from the sun, respectively. Compared with the ejec- 
tion velocities of the 3-4Mq HVSs, those of the 7-15Af0 
HVSs are relatively higher because of their higher mass 
(see Eq. [Ij . Most of the HVSs have proper motions in 
the heliocentric rest frame as large as masyr^^ to a few 
tens masyr"^ (see panel (d) in Figure [S|). These HVSs 
are bright enough to be detected at a distance less than 
a few tens kpc by future telescopes, such as, the Global 
Astrometric Interferometer for Astrophysics spacecraft 
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Fig. 5. — Cumulative distributions of the Galactocentric distance Rqc (panel a), the velocity in the Galactoccntric rest frame (panel 
b), the heliocentric radial velocity in the Galactic rest frame (panel c), and the proper motion in the helio-centric rest frame (panel d), of 
those simulated unbound HVSs (voo > 750km s~^) that were initially associated with the "S-stars". The solid (red), dotted (blue), dashed 
(magenta), and dot-dashed (cyan) curves are for the Unbd-MSO, Disk-MSO, Disk-THO, and Disk-TH2 models, respectively. 



(GAIA), and their proper motions are also large enough 
to be measured. 

Figure |6] shows the sky distribution of the simulated 
HVSs with mass^ 7 — 15Mq, which may represent the 
ejected companions of the "S-stars" , in the Galactic coor- 
dinates by a Hammer- Aitoff projection. In the top panel 
of Figure [HI the positions of the HVSs are projected to 
infinity from the GC, which are consistent with being lo- 
cated close to the CWS disk plane and the Narm plane 
(also projected to infinity) as expected. In the bottom 
panel of Figure [51 the positions of the HVSs are not pro- 
jected to infinity. As seen from the bottom panel, the 
simulated HVSs with mass '--^ 7 — 15Mq lie in the area 
below the projected curves of the CWS disk plane and 
the Narm plane, and most of the se high-mass HVSs re - 
side out of the area surveyed by iBrown et al.l ()2009bD . 
Our calculations show that less than 10% of the simu- 
lated unbound HVSs with mass 7 ^ 15Mq are located 
in the survey area. This may be the reason that none of 
those high-mass HVSs, possibly the companions of the 
"S-stars" , has been discovered in the survey area. If the 
HVSs are initially originated from the CWS disk and 
the Narm plane, the HVSs with high radial velocities are 
also relatively rare in the direction close to the disk nor- 
mals, i.e., {l,b) = (311°, -14°) and (176°, -53°) while 
the HVSs with high proper motions (^ 20masyr~^) is 
relatively numerous in that direction because the veloc- 
ity vectors of HVSs are close to be perpendicular to the 
disk normals. Surveys of HVSs in the southern sky with 
Sky Mapper and others may find such massive HVSs, as 
the possible companions of the "S-stars", and provide 
crucial evidence for whether those S-stars are produced 



by the tidal breakup of stellar binaries. 

Note that the spatial distribution of the HVSs dis- 
cussed in this Section is directly related to the assump- 
tion that the injecting stellar binaries are originated from 
two disk-like stellar structures similar to the CWS disk 
in the GC. However, the other properties of the HVSs 
or the S-stars discussed in this paper are affected by 
whether the injecting stellar binaries are bound to the 
central MBH or not, but not affected by whether they 
are initially on the CWS disk plane or not. 

6. THE INNERMOST CAPTURED STAR 

The Unbd-MSO, Disk-MSO, Disk-THO, Disk-TH2, 
Disk-IMO, Disk-IM2 models roughly produce 147, 111, 
21, 18, 52, and 46 captured stars surviving to the present 
time and with mass in the range of 3 — 7Mq , less massive 
than that of the S-stars, within a distance of ^ 4000AU 
from the MBH (see Tab. [Jj and Section [SJ. The above 
numbers are obtained by calibrating the injection rate 
of the stellar binaries over the past 250 Myr to generate 
17 simulated S-stars similar to the observational number. 
The captured stars with mass in the range of 3 — 7Mq 
could be detected by the next generation telescopes, e.g., 
the Thirty Meter Telescope (TMT) or the European Ex- 
tremely Large Telescope (E-ELT). These low- mass stars 
are potentially important probes for testing the general 
relativistic effects near a MBH, if they are closer to the 
central MBH than S2. In this Section, we estimate the 
probability distribution of the innermost captured low- 
mass stars ('^ 3 — 7Mq) by Monte-Carlo realizations 
based on the calibrated injection rate. 

The left panel of Figure [7] shows the probability distri- 
butions of the semimajor axis of the innermost captured 
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Fig. 6. — Spatial distribution of the simulated unbound HVSs 
with mass 7— ISA/q and Voo > 750km s^^ for the Disk-IMO model. 
The distribution is expressed by a Hammer- AitofT projection in the 
Galactic coordinates. In the top panel, the positions of the HVSs 
are projected to infinity from the GC. The solid red and blue curves 
represent the CWS and NARM disk planes, respectively, which 
are also projected to infinit y. The region above the green curve 
shows the area surveyed bv IBrown et al.l I l2009al ) in the northern 
hemisphere. The open red and blue circles represent the unbound 
HVSs originated from binary stars on the CWS and NARM disk 
planes, respectively. For illustration purpose, all the 58 simulated 
HVSs are shown in the top panel (see Tab.|4]l. The injection rates 
are assumed to be the same for the injections from the CWS and the 
NARM planes and thus the number of HVSs associated with the 
two planes are also the same. In the bottom panel, the positions of 
the HVSs are not projected to infinity. The solid red curves and the 
solid green curve are the same as those for the top panel. The blue 
plus symbols and red crosses represent those unbound HVSs with 
krfl > 275kms~^ injected from the CWS and NARM disk planes, 
respectively. The red open diamonds and the blue open squares 
represent the unbound HVSs with proper motion fi > 5masyr~^ 
injected from the CWS and the NARM disk planes, respectively. 

star with mass ^ 3 — 7Mq resulted from different in- 
jection models. For those models adopting /? = 0, the 
resulted innermost captured star is typically on an or- 
bit with semimajor axis ~ 300AU, and the probability 
that its semimajor axis is less than that of S2 is ^ 99%. 
For the other models adopting (3 = 2, the resulted inner- 
most captured star is on an orbit with semimajor axis 
of ~ 300 — 1500AU and the probability that its semima- 
jor axis is smaller than that of S2 is ^ 60 — 70%. The 
probability to capture a star within the orbit of S2 is 
larger for the (3 = models than for the /3 = 2 mod- 
els. The reason is that relatively more stellar binaries 



can be injected into the immediate vicinity of the central 
MBH and thus more stars can be captured onto orbits 
with smaller semimajor axes (see Tab. |4]) in the models 
adopting (3 = than that adopting (3 = 2. We conclude 
that the probability of a less massive star (3 — 7Mq) ex- 
isting within the S2 orbit is at least 61% and can be up 
to 99%, which may be revealed by future observations 
and then offer important tests to general relativity. 

The right panel of Figure [7] shows the probability distri- 
bution of the pericenter distance of the innermost cap- 
tured low- mass star (~ 3 — 7Mq) resulted from differ- 
ent injection models. For the injection models adopting 
^ = 0, the pericenter distance distribution is concen- 
trated within 50AU; while for the other models adopting 
(3 = 2, the expected pericenter distance is broadly dis- 
tributed over 10 to 200 AU. Nevertheless, the probability 
that the pericenter distance of the innermost captured 
star with mass '--^ 3 — IMq is less than that of S2 (and 
S14) is still significant, i.e., > 55% (or 38%). The in- 
nermost captured star may have its semimajor axis and 
pericenter distance both significantly smaller than those 
of S2, therefore, the GR effects on its orbit may be much 
more significant than that on S2. 

If taking into account the captured stars with even 
lower masses, e.g., IMq, the number of the expected 
captured stars surviving to the present time becomes 
much larger, especially for those models with large 7. 
For example, the numbers of the captured stars with 
mass > IMq are 907, 841, 49, 39, 214, and 117 for the 
six models, respectively. For those lower mass captured 
stars, the semimajor axis and the pericenter of the inner- 
most one could be even closer to the central MBH. 

Some stars may be transported to the vicinity of the 
central MBH by some mechanisms other than the tidal 
breakup of stellar binaries. It is possible that some of 
these stars, with their origins different from the cap- 
tured stars discussed above, exist within the S2 orbit, 
but which is beyond the scope of the study in this paper. 

7. CONCLUSIONS 

In this paper, we investigate the link between the S- 
stars in the GC and the HVSs discovered in the Galactic 
halo under the hypothesis that they are both the prod- 
ucts of the tidal breakup processes of stellar binaries in 
the vicinity of the central MBH. We perform a large num- 
ber of the three-body experiments and the Monte-Carlo 
simulations to realize the tidal breakup processes of stel- 
lar binaries by assuming a continuous binary injection 
rate over the past 250 Myr, and adopting several sets 
of initial settings on the injection of binaries. After the 
tidal breakup of a binary, we follow the dynamical evo- 
lution of the capture d components in the GC by using 
the ARMA model fsee lMadigan et al.l[20To[ ). which takes 
into account both the RR and NR processes, and we also 
trace the kinematic motion of the ejected component in 
the Galactic gravitational potential. 

The properties of the ejected and captured compo- 
nents of the tidally broken-up binaries are naturally 
linked to each other as they are both the products 
of tidal breakup of binaries. For those HVSs discov- 
ered in the Galactic halo with mass '-^ 3 — 4Mq and 
Voo ^ 700 — 1000km s~^, their companions are expected 
to be captured onto orbits with semimajor axis in the 
range ~ 1000 — 8000AU; for the observed S-stars with 
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Fig. 7. — Probability distributions of the semimajor axis [left panel, P(acap)] and the pericenter distance [right panel, P(rp,cap)] of the 
innermost captured star with mass in the range ~ 3 — 7M0 (lower than the masses of the S-stars) at the present time. The solid (red), dotted 
(blue), short-dashed (magenta), dot-dashed (cyan), dot-dot-dot-dashed (green) and long-dashed (yellow) lines show the results obtained 
from the Unbd-MSO, Disk-MSO, Disk-THO, Disk-TH2, Disk-IMO, and Disk-IM2 models, respectively. Note here that these probability 
functions rely on the estimation of the numbers of the simulated detectable captured stars (see Section [Sjl. The estimates are obtained by 
calibrating the injection rate of stellar binaries over the past 250 Myr to generate the same number (17) of simulated S-stars surviving to 
the present time as that of the observed ones. For reference, the position of S2 (or S14) is labeled in the Figure. The probabilities that the 
innermost star is located within the S2 orbit (i.e., Ocap < lOOOAU) are 0.99, 0, 99, 0.99, 0.61, 0.99, and 0.73 for the six models, respectively. 
And the probabilities that the pericenter distance of the inner most star is smaller than that of S2 (S14), i.e., 120 AU (76 AU), are 0.86 
(0.83), 0.81 (0.75), 0.78 (0.70), 0.55 (0.38), 0.80 (0.70), and 0.62 (0.46) for the six models, respectively 



semimajor axis ^ 1000 — 4000AU in the GC, their com- 
panions are expected to be ejected out to the Galactic 
bulge and halo with v^o ^ 500 — 2000km s^^. 

The energy of the captured stars evolves with time be- 
cause of their dynamical interactions with the environ- 
ment. For the captured stars with mass 7 — 15Af0, the 
differences between their present-day energy and their 
initial ones are no more than 30%; for the captured stars 
with mass ~ 3 — 4M0, however, the difference can be 
by order of unity. Therefore, the current semimajor axis 
distribution of the S-stars may provide a good estimation 
on the velocity distribution of their ejected companions 
(e.g., Eq. [5]), but that of the captured stars with mass 
'--^ 3 — 4AfQ does not. The eccentricities of the "S-stars" 
(~ 7 — 15Mo) are close to 1 right after the capture and 
may evolve to low values, and the eccentricity distribu- 
tion of these simulated S-stars at the present time can 
be statistically compatible with the observational ones 
of the S-stars attributed to the RR processes. For those 
captured stars with mass ~ 3 — ^Mq, their eccentricities 
can evolve to even lower values at the present time com- 
pared with the high-mass S-stars (~ 7 — I'qMq) because 
they interact with the environment for a longer time. 

To reproduce both the numbers of the detected HVSs 
and S-stars, the injection rate of binaries need to be on 
the order of 10"^ - IQ-^yv'^ and the IMF of the pri- 
mary components is required to be somewhat top-heavy 
with a slope of ~ 1.6. For the injection models that can 
reproduce the observational results on both the S-stars 
and the HVSs, including the distributions of the semi- 
major axes and eccentricities of the S-stars, the spatial 
and velocity distributions of the detected HVSs, and the 
number ratio of the HVSs to the S-stars, the expected 
number of the ^ 3 — IMq captured companions is ^ 50 
within a distance of ~ 4000AU from the central MBH. 
Future observations on the low-mass captured stars may 
provide a crucial check on whether the S-stars are origi- 
nated from the tidal breakup of stellar binaries. 

The companions of the HVSs, which are captured by 



the central MBH, are usually less massive than that 
of the S-stars (~ 7 — ISMq). The semimajor axis of 
the innermost captured star with mass ~ 3 — 7Mq is 
~' 300 — 1500AU, and the probability that it is smaller 
than that of S2 is - 70% - 90% for the /3 = 2 models and 
~ 99% for the /3 = models. The pericenter distance of 
the innermost captured star with mass ^ 3 — IMq is 
~ 10 — 200 AU and the probability that it is smaller than 
that of S2 (or S14) is also significant, i.e., > 55% (or 
38%). The existence of such a star will provide a probe 
for testing the general relativistic effects in the vicinity 
of a MBH. Future observations by the next generation 
telescopes, such as, TMT or E-ELT, will be able to in- 
vestigate the existence of such a star, and provide im- 
portant constraints on the nature of the central MBH if 
such a star is detected. 

The number of the ejected unbound companions of the 
"S-stars" (see the definition of the "S-stars" at the end of 
Section [5]) is roughly 30 — 40 and the majority of these 
ejected stars are located within a distance of ^ 20 kpc 
from the GC. The number of these ejected companions is 
substantially larger than the number of observed S-stars 
mainly because the observed S-stars are only a fraction 
of the "S-stars" and the rest of the "S-stars" were tidally 
disrupted and do not survive today (see Tabled]). Their 
heliocentric radial velocities in the Galactic rest frame 
range from ^ — 500km s~^ to ~ 1500km s~^ and their 
proper motions in the heliocentric rest frame can be as 
large as ^ 20masyr~^. These high-mass ejected stars are 
bright enough to be detected at a distance less than a few 
ten kpc and their proper motions are also large enough 
to be measured by future telescopes, such as GAIA. The 
majority of the ejected co mpanions of th e S-star s lie out- 
side the area surveyed by iBrown et al.l (|2009af l for our 
observers located at the sun. 
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